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УДК 681.51 
В.Л.ЗАКОВОРОТНЫЙ, ФАМ ДИН ТУНГ 


МОДЕЛИРОВАНИЕ ЭВОЛЮЦИИ ДИНАМИЧЕСКОЙ СИСТЕМЫ, 
ВЗАИМОДЕИСТВУЮЩЕЙ СО СРЕДОИ 


В статье предлагается способ математического моделирования эволюционных из- 
менений в динамической системе, взаимодействующей со средой. При этом пара- 
метры динамической модели среды представляются в виде интегральных операто- 
ров Вольтера относительно траекторий работы и мощности необратимых преобра- 
зований. 

Ключевые слова: моделирование, эволюция, интегральные операторы. 


Введение. Большинство технических систем функционируют в условиях 
взаимодействия отдельных их элементов через различные среды: аэроди- 
намические, гидродинамические, трибологические, технологические и др. 
При этом в пограничных зонах сопряжения некоторых элементов системы 
со средой формируются диссипативные образования, оказывающие влия- 
ние на траектории движений. Свойства этих образований, характеризую- 
щих формируемую динамическую связь, обладают эволюционной изменчи- 
востью. Причём её эволюция зависит от траекторий работы и мощности 
необратимых преобразований, имеющих место в узлах сопряжений элемен- 
тов системы со средой. 

Изучение эволюции технической системы является традиционным, 
и оно связано с проблемой обеспечения надёжности. В частности, при об- 
работке на металлорежущих станках изучение эволюционных преобразо- 
ваний затрагивает такие вопросы, как развитие износа инструмента, изме- 
нение показателей качества изделий в ходе функционирования системы 
резания и др. Однако все эти показатели характеризуют внешнее проявле- 
ние эволюционных изменений. Большее значение имеют эволюционные 
изменения динамической характеристики процесса резания, которая ока- 
зывает влияние на параметры геометрического качества изделия и состоя- 
ние процесса обработки. 
Математическая модель эволюционных преобразований. При по- 
строении динамических моделей используется иерархический принцип, 
основанный на разделении движений на «медленные» движения исполни- 
тельных элементов и «быстрые», которые рассматриваются в вариациях 
относительно траекторий «медленных» движений [1]. Будем считать за- 
данными траектории «медленных» движений исполнительных элементов, 
которые формируются В независимой системе отсчёта 


Х={Х,,Х,,Х,,Х,}' сХ. Здесь и далее символ {---}” означает операцию 
транспонирования. Для определённости, например, применительно к то- 
карной обработке компоненты вектора Х имеют следующий смысл: 
Х | - Координата поперечного перемещения суппорта; Х > - координата 


184 





Вестник ДГТУ, 2006. Т.6. №3(30) 15ВМ 5-7890-0368-0 








перемещения суппорта в направлении скорости резания (очевидно, что в 
традиционной компоновке станка Х, =0); Х. - координата поперечно- 


го перемещения суппорта; Х д - УГЛовая координата положения шпинде- 
ля. Если обработка ведётся с неизменной частотой вращения шпинделя, то 
АХ . /4 = СОП51. Таким образом, в пространстве Х задаются траекто- 
рии исполнительных перемещений станка. Если не принимать во внимание 
погрешности приводов исполнительных перемещений, то эти траектории, 
заданные в виде фазовых траекторий, — есть программа ЧПУ. Например, 
для токарного станка при программировании траектории движения суппор- 
та задаётся скорость подачи на заданном перемещении, то есть фазовая 
траектория {Х з,АХ 3 / Е}. В результате пересечения траекторий про- 
странства Х с расположенной в этом пространстве заготовкой формируют- 
у Т 

ся взаимодействия Р=(Ё (1), Е, (1), Е, (1)) ЕР, заданные своими проекци- 
ями. Эти силы совершают работу при некоторой их мощности. Именно ра- 
бота сил резания при некоторой их мощности приводит к изменению 
свойств процесса резания, вызывая эволюционные преобразования в ди- 
намической системе резания. 

Силы резания зависят от упругих деформаций инструмента относи- 
тельно заготовки. Поэтому введём в рассмотрение векторы упругих дефор- 


маций инструмента Ё= (Е (1), Е (1), Е; ())' е У относительно заготовки 


Е=(Е@), Е, (1), Ё (Е))' ЕЁ. Следовательно, в общем случае силы реза- 


ния являются функциями всех координат пространств Х, \ и 2. Как уже от- 
мечено, в рамках настоящей статьи траектории исполнительных элементов 
станка и геометрия заготовки в пространстве Х считаются заданными. Сле- 
дуя [1], определим уравнение динамики системы, которое дополним эво- 


люционными изменениями параметров р = {Р,, Р2»>° . динамиче- 
ской характеристики процесса резания (в общем случае - среды): 

















2 
г УС От ое. у Е 
ГИ 12 ГИГ ПМ ГИГ Ш 
2 
мн + ср-во, 2.9 ру; 
Ш Ги 11 ГИГ ГИГ Ш 
— (1) 
РА ро-р, | мир (АР) М(РАти =1,2...5; 
0 
1 
АВ = | ма 
0 
№8) =УЕ, (1) 
ле т, М - диагональные матрицы размеромз®з; ве: 


Сс = о — положительно определённые симметричные матрицы 
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жёсткости подсистемы инструмента и заготовки размером 3®3, 
неизменные по координатам перемещения суппорта и неизменные 


при смещении точки равновесия системы; Й = |^, й | Н = |н, в Е 


положительно определённые симметричные матрицы диссипации 
подсистемы инструмента и заготовки размером 3©®3, также неиз- 
менные по координатам перемещения суппорта и неизменные при 


- 1 
смещении точки равновесия системы; У (А-г)=ехр[-—(А- т) 
В 
— ядра интегральных операторов, причём Г» - постоянные работы, 
1 


имеющие размерность км; эти параметры характеризуют эволю- 
ционную наследственность траекторий по мере совершения работы; 


У, — скорость резания, которая в нашем случае считается постоянной. 


Сделаем предварительные замечания по поводу модели (1). В рас- 
сматриваемой системе параметры динамической характеристики процесса 


резания имеют начальные значения Р;0 и значения р; зависящие от 


траектории работы и мощности. В свою очередь, траектория работы и 
мощности является функцией координат состояния системы. В дальнейшем 
будем рассматривать наиболее важный случай, когда работа и мощность 
изменяются настолько медленно, что в пределах одного шага интегрирова- 
ния системы (1) их можно считать неизменными. Другими словами, запаз- 
дывание в вычислении А и М на один шаг не влияет на общую динамику 
системы. Кроме этого в систему (1) входят параметры динамической моде- 
ли подсистем инструмента и заготовки. Методы вычисления этих парамет- 
ров и их идентификации известны [2]. Что касается динамических моделей 
процесса резания, раскрывающих связь координат исполнительных эле- 
ментов станка и упругих деформаций инструмента относительно заготовки, 
то такие модели рассмотрены ранее [2]. В частности, показано, что в 
окрестности стационарной траектории, задаваемой «медленными» движе- 
ниями исполнительных элементов, возможно их линеаризованное пред- 
ставление. Тогда в окрестности стационарной траектории «медленных» 
движений реакция со стороны процесса резания определяется матрицами 
динамической жёсткости и диссипации процесса обработки. Принципиаль- 
ное значение при изучении эволюционной перестройки динамической си- 
стемы имеют параметры ядер интегральных операторов у», (А г), кото- 


рые в данном случае считаются стационарными. Вопросы идентификации 
ядер будут рассмотрены ниже. Вначале, чтобы раскрыть общие принципы 
анализа, будем считать их заданными. 
Особенности эволюционных преобразований динамической си- 
стемы. Рассмотрим систему, обладающую следующими свойствами: 

- силы резания удовлетворяют гипотезе неизменной их ориентации 
в пространстве. Ориентация задаётся угловыми коэффициентами 


ХЕ д.3)! 
- В вариациях относительно стационарной траектории, задаваемой 
Х , динамическую характеристику можно линеаризовать и реакцию со 
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стороны процесса обработки заменить матрицами динамической жёсткости 
и диссипации. Тогда вместо (1) имеем 








у у 
та и +14 и +) = Ел. 
2 
7 
МН + С РИ. 


ЕВЕ (0 +9 (А)(-У @)- 2.0) + ВС“ 1®- о (2) 


А 
Е = 1 [и (А-ОМ(а ат; 
0 


А А 
КА) = оо + ум. (А- маг: ВК) = Ву + п [ур (А- Маг: 
0 0 








1 

Ар) = | мФазм@ =У „ль 
0 

где о. 7, П - согласующие И. имеющие соответственно 


размерность Г 5] эти коэффициенты характеризу- 








ют интенсивность ое а постоянной составляющей си- 
лы резания, жёсткости и диссипации процесса обработки. В этом случае 
эволюционные составляющие постоянной составляющей силы резания, 
коэффициента жёсткости и диссипации соответственно получаются в [кГ], 
[кг/м] и [кГ/м]. 

В системе (2) начальные условия определяются значениями упру- 
гих деформаций инструмента К относительно заготовки 


ь р" на начальной стадии процесса, когда эволюционные 


изменения параметров динамической жёсткости и диссипации отсутствуют, 
то есть @(А)= В(А) =0. Таким образом, 


Уи", © ее 





. 


2(0)={71,25,73}, (3) 





В приведённой модели скорость м. считается величиной постоянной. 
Таким образом, если система имеет точку равновесия асимптотически 
устойчивую, то работу совершает только составляющая силы, направлен- 


ная по скорости резания. Координаты У(0)=Уу” и 7(0) = *, определяю- 
щие начальное состояние системы, являются решениями систем: 


су(0) = х.Е(0);С2(0) = Д.Е (0) =1,23. (4) 


Функции У и ХХ характеризуют смещения точки равновесия, обу- 
словленные особенностями силового взаимодействия на начальном этапе 
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процесса резания и эволюционными преобразованиями системы (2) в связи 
с совершаемой работой при некоторой мощности необратимых преобразо- 
ваний, зависящей от траектории. Старт эволюционных преобразований 
начинается в точке 1 =0, в которой совершённая работа равна нулю. По- 
этому имеет смысл при анализе рассматривать — композицию 


У(А)=У* +УЗ(А) и 7(А)=7* +7С(А), которой во времени соот- 


ветствует композиция У(Р) =У” + ус?) Фи 70=7' +70. Так 
как работа, время, текущие силы и перемещения в (2) взаимосвязаны, то 
можно также рассмотреть функци —У(Х,) Ур ое и 


ИХ) =7”+7С(Х,) смещения точки равновесия по перемещению 


инструмента относительно заготовки, причём здесь Х з - это путь, прой- 


денный инструментом относительно заготовки при изготовлении партии 
изделий. 


Функция А(Г) является возрастающей на всем интервале инте- 


грирования, так как подынтегральное выражение в любой момент времени 
неотрицательная величина. Так как в приведённой упрощённой модели 
матрицы суммарные матрицы жёсткости и диссипации имеют коэффициен- 
ты, которые «медленно» изменяются в ходе эволюционных преобразова- 
ний, то эволюционной является не только точка равновесия, но и изменя- 
ются динамические свойства системы для «быстрых» движений. Поэтому 


эволюционной траектории У“? (г) и 7(г) должны соответствовать тра- 


ектории корней характеристического полинома в комплексной плоскости. 
Однако при высокой интенсивности эволюционных преобразований воз- 
можна такая ситуация, когда сами эволюционные преобразования оказы- 
вают влияние на динамику системы. Интенсивность эволюционных преоб- 
разований в (2) определяется коэффициентами 7 и 7]. 


Приведём пример. Будем считать неизменными динамические ха- 
рактеристики подсистемы инструмента и все условия обработки, кроме 
значений динамической жёсткости и диссипации процесса резания. Основ- 
ные параметры приведены в табл.1, 2, начальные значения жёсткости и 


диссипации с, =100кГ/ мм и ЛО Е. Сила резания на начальном 
мм 


этапе Е, =100к/’; постоянные эволюционной наследственности Г, и 
Тв соответственно равны Т, =50к/м и Тв = 20 кГ’м ; коэффициенты 
ориентации силы: у, = 0,50; у, =0,71; у. =0,50. 

















Таблица 1 
Параметры подсистемы инструмента 
2 кГс 

т кГе .10-3 в ре -103 

2: мм мм мм 
0,25; 0; 0 0,6; 0,1; 0,08 10; 0,3; 0,2 
0; 025; 0 01; 0,5; 0,2 0,3; 1,2; 0,4 
0; 0; 0,25 0,08; 0,2; Ол 0,2; 0,4; 1,6 
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Таблица 2 
Параметры подсистемы заготовки 
2 

й кГс 10-3 АЕ 0 

? мм мм мм 

5,0; 0; 0 5,0; 1,0; 0,8 0,2; 0,1; 0,08 

0; 5,0; 0 1,0; 4,0; 0,6 0,1; 0,2; 0,01 

0; 0; 5,0 0,8; 0,6; 2,0 0,08; 0,01; 1,50 














Мы видим, что в зависимости от параметров интенсивности эволю- 
ЦИОННЫХ преобразований у и п кривые смещения точки равновесия в 


направлении, нормальном к оси вращения заготовки, меняются (рис.1) и 
при некоторых их значениях они напоминают кривые износа инструмента. 
Это лишний раз указывает на связь величины износа инструмента с изме- 
нениями динамической характеристики процесса резания. Заметим, что в 
существующей в настоящее время практике управления точностью обра- 
ботки на основе соответствующих подналадок координат исполнительных 
элементов станков управление осуществляется с помощью непосредствен- 
ного измерения диаметра детали после обработки, что выполнить, как 
правило, сложно, так как требуется дополнительная операция. Приведён- 
ная модель позволяет прогнозировать эти изменения в легко измеримых 
траекториях преобразующей системы станка. 


(79 +2), мм-107' 
14 


1,0 


0,6 


0,2 





Рис.1. Смещение точки равновесия системы в направлении, 
нормальном к оси вращения заготовки, при её эволюционных преобразованиях: 


1-п=0,2, у=01;2-п = 0,2, у=0,5; 3-1=0,2, у=1,0 


Рассмотрим траектории корней характеристического полинома 
(рис.2,а), соответствующие эволюционным траекториям отклонения точки 
равновесия, изображённым на рис.1. Траектории вычислены для характе- 
ристического полинома системы (2) в предположении, что текущие пара- 
метры системы являются замороженными. Рассматриваемой системе отве- 
1 а 
а до 


чают шесть пар комплексно-сопряжённых корней р; и 


р? =-а® = 6 (1=1,2,..-6), расположенных в левой комплексной по- 
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луплоскости, так как замороженная система остаётся устойчивой на всём 
исследуемом отрезке функционирования. Эволюционной траектории каж- 
дого корня или его вещественной и мнимой составляющим однозначно со- 
ответствует эволюционное смещение точки равновесия. На рис.2,6 приве- 
дён пример, иллюстрирующий возможность разделения отдельных эволю- 
ционных стадий на классы по признаку смещения точки равновесия. 


50, с"! 10? 
0,9 


0,8 






И 

ОИ 202-09 ] 
ОИ Е ВЯ 
Гот евочеь 
И |. 


0,5 0,6 0,7 0,8 09 50.с.103 


0,7 











0,6 


0,5 
0,4 








а) 6) 


Рис.2. Эволюционные траектории корней характеристического полинома (а) 
и пример разделения на классы по величине эволюционного смещения 
точки равновесия (6) 


Так как рассматриваемая система является детерминированной и 
эволюционные кривые корней не имеют участков неоднозначности, то 


оценивание УС?) (6) и 72) может быть абсолютно точным. Однако 


возможна постановка вопроса об оценивании эволюционных смещений на 
основе авторегрессионного спектрального анализа с помощью наблюдения 
за вибрационной последовательностью. В этом случае можно сформиро- 
вать информационное пространство, состоящее из вещественных и мнимых 
составляющих корней и использовать методы классификации, подробно 
рассмотренные ранее, например, на основе байесовского классификацион- 
ного правила. 

В примере принята гипотеза о неизменности ориентации сил реза- 
ния в пространстве, которая справедлива лишь в низкочастотной области и 
при малых вариациях износа режущего инструмента. При рассмотрении же 
системы (1) в вариациях, приходится считаться с тем, что в ходе эволюции 
наблюдается перераспределение между составляющими сил, действующих 
в различных направлениях. При этом в суммарных матрицах жёсткости и 
диссипации не только образуются кососимметричные составляющие, но и 
сами матрицы могут стать отрицательно определёнными. 

В этом случае точка равновесия может стать неустойчивой, и тогда 
в пространстве состояния формируются некоторые многообразия [2]. Про- 
иллюстрируем этот эволюционный режим на простейшем примере. 
Будем считать деталь недеформируемой и процесс точения при неизмен- 
ных внешних условиях (скорость, подача, припуск). Для упрощения не бу- 
дем рассматривать влияние работы и мощности на общие силы резания, 
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то есть ЕС?) =0. Таким образом, эволюционно меняются исключительно 
матрицы динамической жёсткости и диссипации процесса. На квазистати- 


ческие позиционные силы в вариациях относительно точки равновесия Х . 
оказывают влияние деформации инструмента в направлении Х! . Эти со- 


ставляющие изменяют суммарную матрицу жёсткости системы. Что касает- 
ся матрицы диссипации процесса резания, то в результате эволюции, как 
уже отмечено ранее, изменяются все её компоненты. Как и ранее, скорость 


резания и две составляющих силы Е ={ В’, Взр будем считать 


постоянными. Таким образом, вместо (3) имеем 


2 
„42, у. и и 

















а? 
Ех 9®,- (ЕС, о, Р(Х. мо, РУХ, ии, 
ах ах ах 
во. О) кучах че ше 3125; 


(5) 


а, (А) =: +а; [ у (А- РМ ат; ВЕСА) = Вог + В ив (А- Маг; 
0 0 


А А 
(= фо: +9, [чи А- Ома ИА) = мод + в [пу (А- ОМемт; 
0 0 


1 
А) = | моасм,(@) =У,|Е> 
0 
Ядра операторов в (5) имеют ту же структуру, что и в (2). Однако 
совершаемая работа соответствует (2) лишь в случае, когда система имеет 
точку равновесия, асимптотически устойчивую. В общем случае в отличие 
от (2) работа и мощность определяются для сил, имеющих различные про- 
екции на оси Х|, Х› и Хз. Выражения для вычисления работы в общем 


случае подробно проанализированы в нашей монографии[1]. Таким обра- 
зом, в ходе эволюции наблюдается перераспределение элементов матрицы 
жёсткости и диссипации системы. При этом в матрицах жёсткости, как пра- 
вило, возрастают те составляющие, которые формируют силовые реакции 
со стороны среды в направлениях Х| и Хз. Коэффициенты матрицы дис- 


сипации варьируются в результате действия следующих факторов: 

- они принципиально зависят от величины запаздывающего аргу- 
мента при формировании изменений сил резания, а запаздывающий аргу- 
мент при неизменной скорости резания увеличивается при возрастании 
объёма пластической деформации в зоне резания, вовлекаемого в пере- 
стройку системы. Поэтому матрицы диссипации при малых вариациях коор- 
динат состояния относительно точки равновесия могут быть отрицатель- 
ными и по модулю они, как правило, возрастают в ходе эволюционных 
преобразований; 

- за счёт кинетической характеристики процесса резания на пада- 
ющем участке зависимости сил от скорости наблюдается эффект отрица- 
тельного трения в области контакта передней грани инструмента со струж- 
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кой и в области, прилегающей к задней поверхности инструмента. Этот 
эффект, как правило, возрастает по мере приработки и развития износа, то 
есть в ходе эволюции системы. 

Рассмотрим пример эволюционного изменения системы резания в 
предположении, что матрицы диссипации остаются неизменными, то есть 
В; =$; =и;=0. Ограничимся случаем, когда работа, совершаемая при 


колебаниях системы в окрестности точки равновесия, пренебрежимо мала, 
по сравнению с работой, совершаемой главной составляющей силы реза- 
ния по мере перемещения инструмента относительно заготовки. Это спра- 
ведливо в тех случаях, когда точка равновесия является асимптотически 
устойчивой или многообразия, формируемые в окрестности этой точки, 
имеют амплитуды, существенно меньшие, чем величина припуска на обра- 
ботку. Тогда вместо (5) необходимо рассмотреть следующую упрощённую 
эволюционную систему резания: 














а2х (и) ах (и _ ах(о_ 
т а + йу ет Е РО }: 
ах@®.__ ах (1) ах (1) Ах т. 
Е(Х, РТ )={Е(Х, р? ), Е>(Х, р? ), Ез(Х, рт )} , 
Е; (Х, = = Ро +а; (А)(-Х1 (1): = 1,2,3; 
А (6) 
а: (=: +а; | и (А- Мат; 


0 


1 

А = | мрак 
0 

№; (@) =У,[Е2 (@) 





) 


(ил Вол); (21 +901); (зд + од) 
где Их =| (№12 + Во.2); (#22 +902); (32 +02) | - суммарная матрица 
(т,з + Роз); (23+ 90,3); (зз + №03) 


диссипации с учётом реакции со стороны процесса резания. 

Рассмотрим конкретный пример для этого случая. Основные пара- 
метры системы приведены в табл. 3. Суммарная матрица диссипации 
условно принята симметричной и положительно определённой. Начальные 
значения сил, определяемые технологическими режимами (1р =2,5мм, 


бр =0,2мм/об, Ур =120,0м/ мин) при точении стали 20Х режущими 
пластинками из Т15Кб равны Юл =40,0кГ, — Ео2 =100,0кГ, 
Юз =60,0кГ'`. Начальные значения коэффициентов жёсткости процесса 
резания соответственно равны «01 =100кГ/мм, 40.2 =80,0кГ/ мм, 
@0,3 = 20,0кГ' / мм. 
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Таблица 3 
Параметры подсистемы инструмента 
с учётом реакции со стороны процесса резания 
2 кг кГ 
т, Е 20— п с. — 10° 

мм мм мм 
0,25; 0; 0 1,2; 0,2; 01 10; 0,5; 0,2 
0; 0,25; 0 0,2; 15; 0,2 0,5; 12; 0,4 
0: 0; 0,25 0,1; 0,2; 14 0,2: 0,4; 16 











Мы видим, что пара корней характеристического полинома после 
некоторой эволюции становятся равными между собой и вещественными 
(точка «А» на рис.3,б). Затем эти корни расходятся по вещественной оси в 
разные стороны, и один из корней пересекает мнимую ось. В этот момент 
«замороженная» система теряет устойчивость точки равновесия, и коорди- 
ната Х,; по закону неустойчивой экспоненты уходит в бесконечность. 


Наблюдается так называемый «подрыв» инструмента. 


ХО, мм 


0,20 


0/15 





0,10 


0,05 








Рис.3. Эволюционные траектории смещения точки равновесия системы 
в направлении Х'\ (а) и корней характеристического полинома (6) для системы (8) 


Система в целом в этом случае претерпевает двойные бифуркаци- 
онные преобразования. Три наложенных друг на друга устойчивых фокуса 
на начальной стадии эволюции преобразуются в точке «А» на рис.3 в два 
устойчивых фокуса, которым соответствуют колебания относительно узла, 
асимптотически стремящегося к точке равновесия. Затем в точке одна из 
траекторий становится неустойчивой, но к ней стягиваются все остальные 
траектории. Приведённый пример показывает, что эволюционные преобра- 
зования в процессе резания характеризуются только износом инструмента 
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и (или) изменением текущих значений диаметра обрабатываемой заготов- 
ки. Это лишь два внешних проявления эволюционных изменений системы. 
Большего внимания, на наш взгляд, заслуживает изменение топологии фа- 
зового пространства, проявляющееся в бифуркационных преобразованиях 
и связанное с ними изменениях динамической характеристики процесса 
резания. 

Идентификация ядер интегральных операторов. При построении 
эволюционных уравнений главную сложность представляет идентификация 
ядер интегральных операторов. Все интегральные операторы в системах 
уравнений являются стационарными. Поэтому возможно их оценивание на 
участке приработки режущего инструмента и установления стационарного 
состояния. Более общий алгоритм идентификации рассмотрен ранее [1]. 
Здесь приведём достаточно простой способ оценивания ядер интегральных 
операторов на примере определения изменений динамической жёсткости 
процесса резания для случая, когда в качестве жёсткости используется 
традиционное скалярное о ней представление, принятое в технологии ма- 
шиностроения. Кроме этого будем считать, что стационарное состояние 
системы является устойчивым. Подчеркнём, что построение эволюционных 
уравнений, позволяющих оценивать смещение точки равновесия, то есть 
изменение диаметра детали, целесообразно выполнять при изготовлении 
изделий, имеющих малую жёсткость, или в тех случаях, когда жёсткость 
подсистемы инструмента мала (например, при растачивании отверстий 
расточными оправками). В этом случае эволюционное уравнение будет 


А 
с) = Хи + [40 + о [зо (А- т)№(г)ат](-У (@) — 21 (1))}; 
0 
А 
Сл = ИЕР +[40 + а [о (А-М()а|(С-" О - 21(1))}; (7) 
0 


т 
А = | Мавме=У ое 
0 
где 0, а; - параметры, подлежащие идентификации; м (А-т)= 





) 


и А-г)} - Ядро оператора, в котором также необходимо 


идентифицировать основной параметр Т.„. 


В процессе исследований заданными параметрами системы (7) яв- 
ляются с1, Су, Х1, 2 И Ур, причём они постоянны. Наблюдаемыми ко- 


ординатами являются У!) (1) и 71(1). Координата У, (г) наблюдаема на ос- 


нове прямого измерения упругих деформаций инструмента, например, с 
помощью тензометрических датчиков. Координата 71(!) определяется на 


основе измерения отклонения диаметра детали после обработки при изме- 
ренной координате У, (г). Измеримой на основе определения упругих де- 


формаций является сила Е в системе, у которой выполняется соотноше- 
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НИе со /с1 =@0/С: =0. Кроме этого по самому смыслу системы (7) при 
А 

#=0 ва [зе (А- маг У (В -71(1)) =0. Поэтому определение &с 
0 

осуществляется на начальной стадии процесса обработки после переходно- 

го процесса, связанного с врезанием инструмента в заготовку, по выраже- 

нию 








С С 
ао = Е У(0)]. 8 
0 ета. 0 я (0)] (8) 


Так как измеримой является У) (1) и Ур =соп5ё , то можно вычис- 
лить фазовую траекторию {А, 4А/ @&= №} и представить её в виде векто- 
ров {А(0), А(1),:--АЮ} и {М№(0), Ма), МЮ}. Кроме этого из (7) спра- 
ведливо 


А 
а, (А) = а [(А-Фм ав, (9) 
0 
О ба, ЗИ. 


У (А)+ 71(А) 2хиУ(А)+ 21 (А)] 

Характерно, что при А=0 ау; =0, так как в этом случае имеем 

2 Ею = < (0) + С. 21 (0) +2 ао [У (0) + 21(0)] - условие статического 
равновесия системы на начальной стадии процесса резания. Поэтому 
смысл функции ау (А) характеризует изменение динамической жёсткости 


процесса резания за счёт совершённой работы необратимых преобразова- 
ний. Так как интегральный оператор является стационарным, и его пара- 


1 

метры ищут в виде м, (А-т)=ехр{ ет. —т)}, то можно сформировать 
[94 

функционал 


[ЕК А АЕ и 
1=Удаз (Ада [Пек ^ Фа? =шт, = (10) 
1=0 0 ы 
позволяющий определить &! и Т„ наилучшим образом в среднеквадрати- 
ческом. Выбор а и Т„, удовлетворяющих функционалу (10), удобно осу- 


ществлять по методу Гаусса - Зайделя. 
Можно ещё в большей степени упростить процедуру оценивания 
@1 и Т» ‚ если ввести в рассмотрение математическое ожидание мощности 


резания с оценкой 


> 


1= 


№=- Ум, (11) 


>| 


1 


|| 
> 


где М = соп51. 
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Тогда функционал (10) существенно упростится: 


А 
1=К ый 

1=У{а5 (Ар - АТ а-е " )}? = ша. (12) 
1=0 


Приведём пример оценивания параметров интегрального оператора 
для случая растачивания отверстия в магний - алюминий - цинковом спла- 
ве МЛ - 5 диаметром 4 = 65,4мм при постоянных режимах резания. Ско- 


рость резания Ур =60,0м/мин, величина подачи на оборот 
$р =0,01мм/об , глубина резания {р = 2,0мм. Инструмент из Т15Кб из 


трёхгранных неперетачиваемых пластинок. Измерение упругих деформа- 
ций расточных оправок в радиальном направлении осуществлялось с по- 
мощью прямого тензометрирования изгибных их деформаций. Пример ха- 
рактеристики функции изменения мощности по работе приведён на рис.4, а, 
а функции изменения ау (.А) - на рис.4,6. По приведённым графикам не- 


сложно — оценить параметры ядра интегрального — оператора. 
Они равны Ть, = 24кГм, @\ =1,08с/кГм?. В этом случае значение эволю- 


ционной составляющей динамической жёсткости процесса резания получа- 
ется в кГ/м. 


М, кГм/с 


755 


250 

















0 50 100 150 АкГи 
а) 6) 


Рис.4. Пример упрощённого способа оценивания параметров ядер интегрального 
оператора: а — фазовая траектория работы-мощности; 6 — функция изменения 
суммарной жёсткости по работе 


При практической реализации систем диагностики и управления 
использование математических моделей для учёта эволюционных преобра- 
зований является недостаточным. Это лишь имитационная модель, позво- 
ляющая определить общую тенденцию изменения траекторий исполни- 
тельных перемещений. Её целесообразно дополнить экспериментально по- 
лученными траекториями эволюции корней характеристического полинома. 

Приведём пример построения таких траекторий для рассмотренного 
случая растачивания отверстий в магний - алюминий - цинковом сплаве МЛ 
— 5 диаметром 4 = 65,4 мм (рис.5). При этом работа сил резания изменя- 


лась от нуля до 300кГм . На диаграммах приведены лишь значимые корни 
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в частотном диапазоне 100 Гу - 2,0кГц. Эти корни формируют простран- 


ство № Й. в котором отображаются отдельные стадии эволюции системы. 
Так как проиллюстрировать графически шестимерное пространство не 
представляется возможным, то на рисунке приведены также простран- 
ственные диаграммы изменения модулей корней р; =а; +Ь; (1=1,2,3), то 


есть |; = \/(а;)? +6." и аргументы, которые в данном случае характе- 


ризуют частоты основных осцилляторов, включённых в АР-модель сигнала 
виброакустической эмиссии. 


О=Ь5кГу 9” 


© = 2,0кГу 








Рис.5. Пример эволюционных диаграмм корней характеристического полинома: 
а- диаграмма эволюции корней характеристического полинома АР-модели 
в комплексной плоскости; 6 — эволюционная диаграмма модулей корней; 
в — эволюционная диаграмма частот 


Ранее показано, что изменения матриц жёсткости и диссипации вы- 
зывают вариации как модулей, так и аргументов. Однако модули, характе- 
ризующие коэффициенты затухания каждого осциллятора, в большей сте- 
пени характеризуют текущее диссипативное влияние процесса резания на 
колебания, а смещения частот (аргументов корней) в основном зависят от 
изменения динамической жёсткости процесса резания. Приведённые здесь 
диаграммы являются типичными для обработки резанием. Особенностью 
рассматриваемых траекторий является достижение вторым корнем единич- 
ной окружности и стягивание на конечной стадии к этому корню третьего 
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(на рисунке показано пунктирной стрелкой). Такое преобразование корней 
обусловлено бифуркациями динамической системы резания после перехода 
второго корня через единичную окружность. В системе в этом случае на 
порядок и более возрастает амплитуда второго осциллятора и подавляются 
колебания третьего. Здесь отметим, что экспериментально полученные 
эволюционные диаграммы характеризуют базу знаний, на которую можно 
опираться для реального оценивания отдельных стадий эволюции. Такое 
оценивание не вызывает сложности, так как опирается на легко измеримые 
временные вибрационные последовательности. 

Выводы. 

1. В ходе функционирования системы резания имеет место её ди- 
намическая перестройка, обусловленная эволюционными изменениями па- 
раметров динамической характеристики процесса обработки. Эволюцион- 
ным траекториям параметров динамической характеристики соответствуют 
траектории корней характеристического полинома в комплексной плоско- 
сти. Наблюдение за траекториями корней в процессе резания удобно вы- 
полнять на основе авторегрессионного спектрального анализа сигнала 
виброакустической эмиссии. Траектории корней, в свою очередь, характе- 
ризуют не используемую ранее информационную базу, позволяющую в ре- 
альном времени диагностировать и управлять процессом резания. 

2. Эволюция параметров динамической характеристики процесса 
резания вызывает изменение параметров многообразий, формируемых в 
окрестностях стационарных траекторий, которые задаются траекториями 
«медленных» движений исполнительных элементов станка. В отдельных 
точках эволюционных траекторий наблюдается изменение топологии фазо- 
вого пространства подсистемы «быстрых» движений, то есть бифуркации, 
которые принципиально меняют динамический режим процесса резания. В 
результате наблюдается смещение траекторий формообразующих движе- 
ний инструмента относительно заготовки, определяющие показатели гео- 
метрического качества изделий. Поэтому наблюдение за траекториями 
корней характеристического полинома позволяет оценивать текущие пока- 
затели геометрического качества изделий непосредственно в ходе обра- 
ботки. 

3. При моделировании эволюционных преобразований в динамиче- 
ской системе резания необходимо использовать функциональные уравне- 
ния, в которых параметры динамической характеристики связаны с траек- 
ториями пространства работа — мощность интегральным преобразованием. 
В результате, с одной стороны, силы необратимых преобразований при ре- 
зании являются зависящими от траекторий, с другой, - сами траектории 
являются зависящими от сил, так как определяются траекториями работы и 
мощности. 
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